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using operators and states 
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We study pure SU(3) gauge theory on a large lattice, using Schrodinger's equation. Our approx- 
imate solution uses a basis of roughly 1000 states. Gauge invariance is recovered when the color 
content of the ground state is extrapolated to zero. We are able to identify the gauge invariant 
excitations that remain when the extrapolation is performed. In the weak coupling limit, we obtain 
promising results when we compare the excitation energies (masses) to known exact results, which 
we derive. We discuss the application of our nonperturbative method in the regime where glueballs 
are present. 

PACS numbers: 11.15.Ha, 12.38.Gc, 11.15.Tk 



I. INTRODUCTION. 

Pure SU(3) gauge theory on the lattice can be for- 
mulated in Hamiltonian form. Matrix elements of the 
Hamiltonian using any finite orthonormal basis can be 
used to form a Hamiltonian matrix, and the eigenvalues 
of this matrix, ordered in increasing energy, are each an 
upper limit to the corresponding eigenvalue of the origi- 
nal Hamiltonian. If the basis is provided with parame- 
ters, these can be used to lower the approximate energies 
closer to their true values. 

Variational calculations run into problems unless the 
basis states are carefully chosen. This is particularly true 
when many degrees of freedom (DF) are involved. The 
most serious problem is that the eigenvalues of the matrix 
may be useless because they lie far above the eigenval- 
ues of the Hamiltonian and do not accurately portray the 
physical implications of the Hamiltonian, quantitatively 
or even qualitatively. Another problem associated with 
many DF is that the basis size can become unmanage- 
able. If we work on a L x L x L spatial lattice, there 
are 24L 3 degrees of freedom. (There are 3 links emerg- 
ing from each site on the lattice, and 8 colors per link.) 
When we introduce two states for each DF, as we will, we 
have a basis of 2 24i states. Further cuts must be made. 
As we do this, we also want to build symmetries of the 
Hamiltonian into the basis so that quantum numbers can 
be assigned to the eigenvalues of the matrix. The founda- 
tion of a variational calculation is therefore the crafting 
of a well-chosen basis. 

The main purpose of this paper is to construct a basis 
for pure SU(3) lattice gauge theory and demonstrate that 
it is well-chosen. In this Introduction we outline the is- 
sues that must be faced in the construction process, and 
how we cope with them. In following sections we com- 
pare results obtained with our basis to known results, 



thus testing the adequacy of the basis. This compari- 
son can be made only in a limiting case, but that case 
imposes significant tests. The methods we use can be ap- 
plied away from the limit to the "physical" case, where 
the lattice is inhabited by glueballs. We outline how to 
do this at the end of the paper. 

The approach we take in this paper is clearly an ap- 
proximation. On the other hand, the customary evalua- 
tion of the path integral by Monte Carlo sampling can, 
in principle, be arbitrarily accurate as more computing 
resources are devoted to the project. But that project 
has been ongoing for 30 years, using the most advanced 
computing facilities in the world. The payoff to our ap- 
proach is that the computer resources required are very 
modest. 



A. The Hamiltonian 

We begin with a description of the elements required 
to write the Hamiltonian for pure SU(3) local gauge the- 
ory on a three dimensional lattice. We label links on our 
lattice by the pair (s, m), where s is the site from which a 
link departs, with < s x , s y , s z < L = 0, and the lattice 
is an L x L x L cube, m = x, y or z is the direction of 
the link. The variables assigned to a link degree of free- 
dom specify a point on the SU(3) manifold, a compact 
manifold of eight dimensions labelled by a = 1,...,8, 
the color index. In the Hamiltonian we use the funda- 
mental representation matrix U StTn , which depends on the 
point on SU(3) for link (s,m). We also use the opera- 
tors Ji,a(s,m) with a = 1,...,8. The action of Jh,a 
on a representation of SU(3) is to left-multiply it by the 
ath generator matrix for that representatrion. Similarly, 
J Ria acts to right-multiply. J 2 = J2^i, a = J2^R, a 
multiplies by the quadratic Casimir operator for SU(3). 
The Hamiltonian is then: 



'Electronic address: bronzan@physics.rutgers.edu 



H A = ^ff(s), 

S 



2 



(1) 



x J2 [ Tr (u s , m u s+ih , n ul +ilM ul n ) + H.c. - 6] . (2) 

m<n 

-Ha is the Kogut-Susskind Hamiltonian. [2J It commutes 
with the generators of local gauge transformations at site 
s: 

Qa(s) = -JL, a {s,X) - JL,a(s,V) ~ Jl,o,{ S > Z ) 

+JrA s -*,«)+ J R , a {s - y, y) 

+jR,a{s - Z,z). (3) 

The Hamiltonian formulation of lattice gauge theory 
has been used for a number of studies, mostly based on 
loop expansions that differ from the approach we take 
here. For a recent review, see 0. 



B. States for Each Link. 

The first step in building a basis is to choose states on 
the SU(3) manifold for each link. To make the choice, we 
need to know where on the manifold the eigenfunctions 
of Ha have support (are large), and this depends on the 
coupling g. We begin with small g. It is convenient to use 
the parametrization for the fundamental representation 
matrix 



U = aip(ig\ a R a /2) 



(4) 



where A a is a Gell-Mann matrix. |4j We next develop as- 
cending power series in g for the elements in the Hamil- 
tonian 



Id i d 
jL ' a ~ ~gdR a + 2 UbcRb dR c + ' 

1 d _ * d 
jR ' a ~ igdR a ~ 2 fabcRh dR c + ' 

2 

U = l + |(A Q i? a )-^(A a i? Q ) 2 



(5) 



To study support, we examine a "lattice" consisting of 
one link. To lowest order in <?, the Hamiltonian is 



R 2 



*-T^-?P™-«*--5s!f + -5 



(6) 



If we treat R as an eight dimensional Euclidean vector, 
the ground state of H' is 



t/> = cxp(-R 2 /2). 



(7) 



The Euclidean approximation is justified when g is small 
because then the group parameters a a — gR a are small 
where the wave function is large, and the fact that SU(3) 
is not Euclidean can be ignored. At weak coupling, there- 
fore, wave functions have support only near the group 



origin. As g grows, the region on the manifold where the 
wave function is large expands. At large g we know from 
strong coupling expansions that wave functions have sup- 
port everywhere on the manifoldf|. Therefore, we seek 
Gaussian-like states on the manifold, centered on the ori- 
gin, with width that can be tuned as a variational param- 
eter. 

Gaussian-like states on the group are generated by the 
kernel of the heat equation^] . 



dt 



J 2 



(8) 



The weak coupling, Euclidean approximation demon- 
states that the kernel is Gaussian: 



dt 



1 

7 



(R 1 |e iv2 ^ 2 |R 2 > = - 



- 9 2 (R!-R 2 ) 2 /4t 

{^t/giy 



(9) 

We will use the scaled width variable r = t/g 2 . Com- 
paring Eqs. Q and iJSJ, for one link at weak coupling 
t = 1/2. This is the order of magnitude for r that we 
will find for a large lattice at weak coupling. 
The heat kernel on SU(3) is0 

i/> ai (a) = ({a}\e- tj2 \{ ai }) (10) 

= Y, d ^l)^ tHP ' q) X (P - q HM{a 1 }- 1 ), 
p,q 

where {a}, {ai} are group parameters, the char- 

acter, d(p, q) is the dimension of representation (p, q), and 
A(p, q) is the Casimir eigenvalue, [a, Note that when t 
is small (narrow Gaussian), many representations con- 
tribute. This mirrors the situation in Euclidean space, 
where many harmonics are required to build up a narrow 
Gaussian. 

In Eucliean space, Gaussians can be extended to a com- 
plete set of basis states by multiplication by polynomials 
or by differentiation, and the same can be done on SU(3). 
The simplest of these states are 



4> = lfjal(a)\a 1= 0, 
4>a = JLa(ai)lp ai ( a )\ai=0- 



(11) 



The orthogonality properties of these states areQ 

{4>\4>) = Si; (0 O |0) = 0; (<t>a\<l>b) = S ab S 2 ; 
5i = ]Td 2 (p, 9 )e- 2tA ^; (12) 



S 2 



p.q 
1 



^d 2 (p,q)\(p,q)e- 2tx ^. 



Matrix elements of the operators in Eq. among states 
<j) and 4> a are computed in Ref. 6]. The dynamics of 
SU(3) is invoked through the use of these matrix ele- 
ments. 

On each link, we allow single excitations to occur in 
each of the eight colors. We exclude states that are mul- 
tiply excited in the same color direction. Since each of 
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the eight color directions on a link can be unexcited (4>- 
like) or singly excited (</> a -like), our basis has 2 8 = 256 
states per link. An important topic of the paper will be 
to test the adequacy of limiting the basis to 256 states 
on each link in this way. 

States in which two color directions on a link are ex- 
cited, designated (f> a b, are studied in Ref. |6J, but states in 
which 3, . . . , 8 directions are excited have not been stud- 
ied. Extending the number of color directions excited on 
a link in the explicit manner of Ref. is laborious and 
will not be attempted in this paper. The process is also 
ambiguous because the operator Jha used in the defini- 
tion of 4> a and 4> a b can be replaced by others, like J7r q - 
Instead, in this paper we introduce states with "multi- 
color" excitations on a link, and matrix elements involv- 
ing these states, in a straighforward way we describe in 
the next subsection. 



C. States on the Lattice 

The lattice basis constructed by putting each link in 
one of the 256 states we have introduced has (256) 3L = 
2 24i elements, which is far too many to handle. In addi- 
tion, local degrees of freedom constitute an awkward ba- 
sis for describing phenomena that range over many links. 
Continuum field theory instructs us how to depart from 
local degrees of freedom by recasting the theory in mo- 
mentum space. We adopt that method in this paper. 

We introduce operators A a (s,m). These are bosonic, 
satisfying the usual bosonic commutation rules. The 
ground state |0) corresponds to the normalized state 
Ify/y/Si on a link, and —iA a \0) corresponds to the state 
\<t>a) /V^2 on the link. At this point we can describe how 
we represent a normalized state in which two distinct col- 
ors are excited on the same link; it is (— iA a )(— iAUO), 
where a ^= b. 

We now follow field theory and introduce nonlocal 
states by passing to momentum space, introducing the 
bosonic operators 

S (k, m) = ^^ e - 2 " k -( s +*/ 2 )/M a (s, m). (13) 



D. The Basis Generator G 

We return to the requirement that the state on each 
link be required to be in one of 256 states. Applications 
of operators (— iA\) k with k > 2 produce states beyond 
the ones we've introduced. We refer to such states as 
"overexcited" . Overexcited states will appear unless we 
suppress them. First we introduce an operator that mea- 
sures the excitation of the color/link degrees of freedom. 
It is 

N= Y, Al(s,m)A a (s,m), (14) 

s,m,a 



A color/link degree of freedom in the k-th bosonic state 
is an eigenfunction of N with eigenvalue k. Our method 
of minimizing the contribution of overexcited states to 
our basis is to control the expectation value of N in our 
basis states. 

With the introduction of N, we can begin to construct 
the basis. First we write H A in terms of our bosonic op- 
erators. Coefficients are chosen so that matrix elements 
on the basis of the two states for each color/link agree 
with matrix elements derived in Ref. Next we de- 
compose H A into a term H A0 , quadratic in A and A' 
and a remainder: 

H A = H A . + H AA (15) 

Equipped with this decomposition, we construct a gen- 
erator operator G, whose low-lying eigenstates will com- 
pose our basis. The first expression for G is 

G A = H A>0 + £A. (16) 

By limiting G A to quadratic terms, it becomes an op- 
erator that can be diagonalizcd. Its eigenstates can be 
constructed, leading to an explicit basis. We take £ to be 
a positive parameter, so the term in G A reduces the 
contribution of overexcited states to the low energy spec- 
trum of G A , where the states in our basis are found. By 
changing £, we can adjust (N) for the eigenstates of G A . 
Define N aV e to be the expectation of N per color/link DF 
in the ground state of G A : N ave = (0| 7V|0) /(24i 3 ). Most 
of our simulations will be carried out with N ave = 1/2. 
This small value suggests that the probability of find- 
ing a color/link to be overexcited is small, and it is. 
We will present a formula giving the probability that a 
color/link is in any overexcited state, and the typical re- 
sult (when Nave = 1/2) is about 1/8. Therefore, overex- 
cited states contribute to results computed in our basis, 
but not much. We have some ability to check on the 
influence of overexcited states by changing £ and hence 
Nave- We do not want to make N ave too small, however, 
because we do not want to suppress the 256 states we 
allow for each link. Another method of excluding overex- 
cited states would be to map pure SU(3) onto a spin 
model, with each color/link being in either a "spin up" 
or "spin down" state. This would eliminate overexcited 
states altogether, but it would deny us the methods of 
field theory, particularly the simple introduction of non- 
local degrees of freedom by transforming to momentum 
space. 

To summarize our procedure, we have adopted a 
bosonic formalism. This allows us to pass to momentum 
space and construct a basis that is nonlocal. The price 
we pay is that there are contributions from overexcited 
bosonic states that are unwanted. By taking our basis 
states to be eigenstates of the generator G A , we can use 
the parameter £ to minimize and assess the effects of the 
overexcited states. 
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E. Gauge Invar iance 

We seek states having Q a (s)\il>) = 0. There are three 
options for excluding non-gauge invariant states. One 
is to rewrite the theory in terms of gauge invariant DF 
only. Unfortunately, the resulting Hamiltonians are non- 
local and are not manifestly under rotations of the cu- 
bic lattice 0, They are unsuited for approximation 
schemes. The second option is to use Hamiltonian (QJ 
but allow only gauge invariant states in the basis. This 
is done when loop expansions are employed. [3j] However, 
sound arguments have led us adopt the eigenstates of Ga 
as basis states, and these states are not gauge invariant. 

The third option is to compute with a nonzero expec- 
tation value, per site, of the group quadratic Casimir 
operator: Q 2 ave = (0|Q 2 |0)/i 3 . Here |0) is the ground 
state of Ga, and 

s . a 

This is the option we adopt. We will calculate results 
for several values of Q 2 ve and extrapolate to the limit 
Qave = at the end. This was done successfully for the 
ground state in Ref. and we use the same method 
here. Note that Q a (s) 2 is a positive operator, so when 
Qave = 0; g au g e invariance is imposed at every site on 
the lattice. 

Once we include gauge non invariant states in our ba- 
sis, we want to control their number. This can be ac- 
complished by adding a term to G that raises the eigen- 
value of such states. To do this, we decompose Q 2 into 
quadratic and higher terms, just as we did for Ha 

Q 2 = Ql + Q\- (18) 

We add term proportional to Qq to our generator, which 
we call Gb after the modification. 

G B =G A + V9 2 Ql = H A .o + + V g 2 Q 2 . (19) 

The factor g 2 is included because g 2 Q 2 ve is finite as g — > 
0. 

Finally we introduce Hb, the Hamiltonian we use in 
the rest of this paper. 

H B =H A + V9 2 Q 2 - (20) 

Ha and Hb have the same set of physical, gauge invari- 
ant eigenstates. However, the non-gauge invariant eigen- 
states of Hb are raised in energy by the positive operator 
i]g 2 Q 2 . By varying r\ we can control the number of gauge 
non invariant states in the low energy spectrum. 

Our basis states are the low-lying eigenstates of gener- 
ator Gb- To diagonalize Gb we must remove operators of 
the form B^B^ and BB. This is accomplished by means 
of the Bogoliubov transformation to new bosonic opera- 
tors E: 

B«(k,m)= (21) 
cosh0(k).E o (k,m) - sinh 0(k)£^ (L -k,m), 



with 0(k) = 0(L-k) and L = (L, L, L). (Directional 
indices for will be introduced later, after we transform 
to the polarization basis.) We use the vacuum state of 
Gb after the Bogoliubov transformation to determine the 
expectation values N ave and Q 2 ve - These parameters are 
set primarily by choice of £ and r. r\ is used to set the 
number of non gauge invariant states in the low energy 
spectrum of Gb. 

F. Group Theory, Degeneracy and Spin 

Hamiltonian is invariant under the group of rota- 
tions of the cubic lattice. Associated with this symmetry 
is degeneracy; states belonging to different rows of the 
same irreducible representation of the group are degen- 
erate. Such degeneracy complicates computer work. It 
can be reduced by forming linear combinations of our ba- 
sis states that belong to the same row of an irreducible 
representation of the group: 

m=^V ( i ? I )*(R)U(R)\a), (22) 

R 

where is an irreducible representation of the group, 
and U(R)\a) is the image of \a) under a rotation of the 
cube. 

The group of rotations of the cube has 24 elements[To|. 
It has five irreducible representations (because there are 
five classes), which we label by their dimensions: 1a, 
1b, 2, 3a and 3b- (The sum of the squares of these 
dimensions is the order of the group, implying that Eq. 
(|22|1 provides 24 states for each \a).) 

An additional reason to take account of the symme- 
try group is that we wish to identify the spins of the 
states we find. Since our lattice breaks the 0(3) sym- 
metry of continuum field theory, this identification is not 
obvious. Note, however that each of the 24 rotations of 
the cube corresponds to an element of 0(3). The iden- 
tity element corresponds to a rotation by 0°, and each of 
the other elements corresponds to a rotation by 90°, 120° 
or a multiple of these about an axis through the center 
of the cube. It follows that every irreducible represen- 
tation of 0(3) generates a representation of the group 
of the cube, which is generally reducible. The first few 
Clebsch-Gordon series giving the content of 0(3)-induced 
representations are 

(spinO) — ► 1a, 

(spinl) — > 3 A , (23) 

(spin2) — ► 2 3 B , 

(spin3) — ► 1 b (B3a®3 B - 

In these four cases, a representation of the group of the 
cube is identified with a spin. However, an eigenstate be- 
longing to the 2 representation of the cube can be associ- 
ated with spin 2 only if there is an eigenstate belonging 
to the 3b of nearly the same energy. If this is not the 
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case, it may not be possible to identify the 2 eigenstate 
with a continuum state. 

In this paper we will study spin states. T>(R) = 1 for 
this case, and 

\^)=Aj2u(R)\a) (24) 



II. WEAK COUPLING 

We need a known reference theory which we can use 
to evaluate the adequacy of the approximations laid out 
in Sec.[l] For this pupose we use the weak coupling limit 
of the theory. This is the limit in which the correla- 
tion length has grown beyond the lattice length L. This 
regime is reached for couplings g below those in the scal- 
ing region in which glueballs are dynamically contained 
within the lattice. The weak coupling limit is unphysical, 
but it provides a nontrivail test of the numerical accuracy 
of calculations using the basis produced by state genera- 
tor Gb- In the weak coupling limit we can also test how 
well we can spot gauge-invariant excitations in a con- 
geries of states that are excitations of Q 2 . Finally, we 
can test the extrapolation to Q 2 ve = 0. 

To obtain the weak coupling form of the theory, use the 
weak coupling operators of Eq. (JSJ) in the Hamiltonian, 
Eq. QJ. We retain only the leading terms in powers of 
g, which produces a free field theory: 



We find that 



1 y 

2 ^ dR a (s,m) 2 

s,m,a ■ ' 

+i ^2 [R a (s,m) + R a (s + m,n) 

s,m<n,a 

-R a (s + n, m) - R a (s, n)] 2 . 
The generators of local gauge transformations are 



(25) 



?W,o(s) = — 



d 



d 



dR a (s,x) dR a (s,y) 
d d 



dR a (s,z) dR a (s-x,x) 
d 8 



8R a (s-y,y) dR a (s-z,; 



(26) 



We diagonalize this theory by introducing bosonic op- 
erators. 

R a (s,m) = y /a[A a (a,m)+Al(s,m)], (27) 
d 1 



dR a (s,m) 



[A a (s,m)- Al(s,m)]. 



a is an arbitrary parameter, and physical consequences of 
Hamiltonian 125|) do not depend on a. Next we introduce 
bosonic operators in momentum space: 

A a (s, m) = ^J2 e 27 " k - (s+ * /2)/i 5a(k, m) (28) 



R a (s, m) + R a (s + m, n) - R a (s + n, to) - R a (s, n) 
= 2^^ e27rik ' S/Le " (fcm+fc " )/L ( 29 ) 



su , [ ?^HL J W a (k,n) -sin (^j^ ) \Y„\k. m i 



where W a (k,m) = B a (k, to) — B\(L — k, to). This may 
be written concisely in terms of the unit longitudinal po- 
larization vector 



e (3) (k) =w(k) 
irk z 

T 



Trk x 



ysm 



(30) 



: (k) = £sin 2 



and the vector 

W„(k) = xW„(k, x) + yW a (k, y) + zW a {k, z). (31) 

as 

R a (s, to) + R a (s + m, n) - R a (s + n, to) - R a (s, n) 



X UJ 



(k) [f< 3 )xW a (k) 



(32) 



with m, n, p a permutation of x, y, z. The sign of the 
term is positive (negative) for an even (odd) permutation. 
Introduce the real transverse polarization vectors e*- 1 ^ (k) 
and e 2 )(k) so that the three polarization vectors form 
an orthonormal right-handed basis. The vectors are cho- 
sen so that e( s )(k) = e^ s ^(L — k). The polarization basis 
degrees of freedom are 



BW(k)=£eW(k)B Q (k,TO). 



(33) 



Then 



2 X] [Ra(s,m) + R a (s + m,n) 

s,m<n,a 

— i?(S + n, to) — R a (s, n)] 2 

2 

= 2a^^^ 2 (k)^)(k)^)t(k). 

k.a s=l 

Altogether, the weak coupling Hamiltonian is 



(34) 



k,a s— 1 


hfl«t(L 


-k)] 


x[i?i s)t (k) 


+ BW(L 


-k)] 


+2a]T^> 2 (k)[i^)(k)- 

k,a s=l 


-B«t(L 


-k)] 


x[Bi s)t (k) 




-k)] 



(35) 
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We now use the Bogoliubov transformation, Eq. (|21(l . 
on the transverse modes, s = 1,2. The parameter is 
29^ (k) = - ln[4<7w(k)]. The "transverse" mode Hamil- 
tonian becomes 

2 

H W . T = L 3 £ + 2^^c(k)^)t(k)i?( s )(k), (36) 

k,a 2—1 

and the ground state energy density is 



£b = 55>(k). 



£o depends weakly on the lattice size. When L = 11, 
£ Q = 19.0999, while when L = 31, £ = 19.1008. 

The longitudinal modes (s = 3) cannot be brought 
into this form. A glance at Eq. 1)25(1 shows that the 
free field theory is a set of coupled oscillators. How- 
ever the "restoring" term vanishes for the longitudinal 
modes, which therefore have a continuous spectrum of 
eigenstates corresponding to free motion. The longitudi- 
nal mode Hamiltonian is 

^ = ^E[^ 3) ( k )+^ 3)t ( L -^ 

[B^(k) + B^(L-k)] . (38) 



k.a 

X 



Define four hermitian operators associated with mode k 
in this sum: 



Pi,. = h-B^M+B^M 

At 

-fl( 3 )(L-k)+£< 3 >t(L-k)], 
P2 , a = i[s( 3 )(k) + B( 3 )t(k) 

-B( 3 )(L-k)-B( 3 )t(L-k)], 
= i[si 3 )(k) + B( 3 )t(k) 

+5( 3 )(L-k) +J B( 3 )t(L-k)], 
n 2 , a = i[5( 3 )(k)-B( 3 >t(k) 

-S( 3 )(L-k) +J B( 3 )t(L-k)]. 
These operators satisfy canonical commutation rules. 
[Pi.a, P2,b] = 0, [7r ljQ , ir 2 , b ] = 0, 



(39) 



(40) 



In terms of these operators, the longitudinal Hamiltonian 
is 

H ^ = r- E^ 3) ( k ) + B ^( L - k )] 



k , a 



x[5( 3 )t(k)+i?( 3 >(L-k)] (41) 
= ^E{[^( k )] 2 + [^a(k)] 2 }- 



Thus we need eigenstates of 7Ti jQ and ii2,a- Since 
[fimij,] = 0, we can look for common eigenfunctions: 



7Tl, a |pi,p 2 ) =Pl|pi,P2>, 7T2, a |pi,P2) =P2|pi,P2> 



In Fock space the eigenfunction is 



(42) 



\P1,P2)= £ ^[^ 3)t (k)] m [si 3)t (L-k)]™|0), 

^ ' 777. ' 77 ' 



m.n—0 



(37) with 



(0 



777+77 



(43) 



dAe-^ Api m (A+p 2 /V2) (44) 



x^„(A- P2/V2). 

The functions in this expression belong to a set of or- 
thonormal harmonic oscillator eigenfunctions: 



</>m(A) 



1 







1 



A 



-X 2 /2 



(45) 



The eigenvalues are continuous, with —00 < P\,p% < 00, 
and the states are orthonormal: 



(puPzlPiifc) = $(pi -h)HP2 -P2) 



(46) 



The mean excitation is N ave = 00 for the weak cou- 
pling theory. (A simple way to see this is provided by 
Eq. I|70l) . First set 77 = and then examine the limit as 
£ — * 0.) This means that overexcited states contribute 
significantly to the eigenstates. This result highlights 
an important difference between the ground state of the 
weak coupling theory and the ground state of generator 
Gb, where we set N ave = 0.5 in most of our simula- 
tions. Our approximation scheme is viable only if we can 
demonstrate the accuracy of calculations using our basis 
of states, despite this difference in N ave . 

The operator Q 2 may also be written in terms of the 
longitudinal operators we have introduced. 



Q 2 = ^E- 2 ( k )[^ 3) (k)+Si 3)t (L-k)] 



go- 



k.a 



;[Bi 3 )t(k)+i?i 3 )(L-k)] 
-— ^> 2 (k){[7r 1)£l (k)] 2 



(47) 



q 2 cr ' 



k2,a(k)] 2 }- 



Since uj 2 < 3, we see that the ground state energy 
density for general Q 2 ve is 



24 



(48) 



k.a 



Energy levels above the ground state are comprised of a 
continuous spectrum of excitations of Q 2 and a discrete 
spectrum of gauge-invariant excitations whose energies 
are sums of 2a; (k) for different values of k. 
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Pure SU(3) is believed to exhibit color confinement, 
but our free field theory does not, nor will the theories 
we propose for g 2 > 0. We contend that this is accept- 
able as long as we compute low energy properties of the 
theory such as glueball masses. (Computations of the 
spectrum of charmonium show that the confinement po- 
tential plays a minor role there.) However, we must take 
account of color confinement by always contracting color 
indices a into guage-invariant combinations. For glue- 
balls at rest, we will use a basis of gluon pair states of 
the form J2 a E { a s)ji (k)E d s)1< (L - k)|0). (In this paper we 
replace negative wave vectors — k by L — k so that all 
vector components lie in the interval < < L.) It 
will be seen that such states arise naturally in the course 
of calculation. Thus, low-lying gauge invariant states at 
rest have energies (masses) 4u(k). This is the result we 
hope to replicate using our approximate eigenstates of 

The spectrum of transverse excitations is highly de- 
generate. Some of these degeneracies depend on familiar 
trigonometric identities and the trigonometric form of 
w(k). For example, when L is even, and k y — L/2 ± k x , 
uj 2 = 1 + sin 2 (nk z /L) for any k x . When L is divisible by 



3, and we take k v 



then 



(L/3) - k x , and k z = (L/3) + k x 



1.5 for any k x . These degeneracies are lattice 



artifacts. We conjecture that they are not present when 
L is a prime number, and we have verified that they are 
absent when L= 11, 31 or 53. When removing degener- 
acy is an important consideration, we will choose L to be 
a prime number. 

Further degeneracies are associated with symmetries of 
the cubic lattice. The factor Lu(k x , k y , k z ) is unchanged 
when the momentum components are permuted (factor 
of 1, 3 or 6) or k^ is replaced by L — k^ (factor of 8). 
There are two transverse polarizations, so up to 96 tran- 
verse polarization states have the same energy. The de- 
generacy of transverse pair states is up to 48. When 
states are grouped to form representations of the rota- 
tion group of the cube, these states are organized into 24 
orthogonal linear combinations of differing "spin" . This 
leaves a two-fold degeneracy, and we expect the trans- 
verse polarization spin zero states we study to be doubly 
degenerate. 

We can choose polarization vectors so that the remain- 
ing double degeneracy is labeled by the polarization in- 
dex s = 1,2. We generate states that are representations 
of the rotation group of the cube by specifying a "seed" 
wave vector in the interval 

< k x < k y < k z < [L- l)/2; (exclude k = 0.) (49) 

Then the sets {Rk} and {L— Kk} span all L 3 — 1 wave 
vectors as R ranges over the 24 elements of the rotation 
group and k ranges over the set of seed vectors. For 
vectors in Eq. I|49|) we choose ejj (k) so that e x 1 ' (k) = 
and ^(k) > 0. = e'^xe' 1 *. Then choose 



These choices assure that all states contributing to a rep- 
resentation of the rotation group have the same polariza- 
tion quantum number. The index s becomes a quan- 
tum number that distinguishes the two degenerate spin 
states at each energy eigenvalue. 

The results we expect for pure SU(3) differ from these. 
The operators Q a (s) satisfy the Lie algebra of SU(3), so 
the spectrum of Q 2 is discrete, consisting of sums of in- 
teger multiples of the quadratic Casimir eigenvalues of 
SU(3): \{p,q) = 4/3, 8/3, .... Thus g 2 Q 2 , which ap- 
pears here, has a discrete spectrum, but one whose spac- 
ing goes to zero at g = 0. The disappearance of the 
gaps between states rationalizes the continuous spectrum 
we have found for g 2 Q 2 when we use the weak coupling 
Hamiltonian. Note that in the weak coupling limit the 
eight gauge generators at each site commute, so the lo- 
cal gauge invariance group becomes Abelian, and is no 
longer SU(3). 



III. SU(3) HAMILTONIAN USING OUR BASIS 

The basis states we use are Fock states generated by 
applying products of operators E„ (k) to the ground 
state of generator Gb, Eq. (|19fl . We write Hamiltonian 
Ha, Eq. Q), as an operator on this basis, using the 
matrix elements of Ref. @ . The correspondence between 
the Fock states and the states in Rcf. 6] is 

|0>~|O)v^, |0 o >~-i4|O}^. (51) 

For example, the four matrix elements among |0) and 
\4> a ) are correctly given by the operators 



Jh(R),a — — 




4)+l(-%./W^K, (52) 



where repeated color indices are summed. 

In the Hamiltonian we also need the operator J 2 . Ma- 
trix elements of this operator are given in Ref. 0, but 
we can also use the relations J 2 = Jh.aJh.a = Jn,aJR,a 
and Eq. I|52|) to compute the operator. The expressions 
differ because when we begin with Eq. (|52|) , intermediate 
states between the two factors do not form a complete set. 
A choice must be made for J 2 , and we will use the square 
of the operator in Eq. 1)521) in this paper because that is 
what appears naturally in the operator Q 2 — Q a Q a - 

J" = — (A a — ^.q) + -A a A a — -fpqafrsaA^AlAqAs. 

The operator reproducing the four matrix elements of 
U among and \4> a ) is 



U = 



Si A^B~ 2 



Xc(A c + Al) 



(54) 



e (s) (i?k) = eW(L-iZk) = i?e ( ^(k). 



(50) 



+ 



DKo 
2S 2 



Ko 
Si 



AtA c + -^{X c ,X d }AlA d 
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The matrix element of the operator W on our basis is 
obtained by changing the sign of the second term in this 
equation. Expressions for the matrix element coefficients 
from Ref. [6j are: 

1 



Kn = 



^e(p, q)f(p, q), 



3 ^ 



K, = - 



4 ^ 



DK Q 



3^ 

p,q 



e(p,q)f(p, q) 



e{p, q)f(p, q) 



i(p,q) . i(q,p) 



20 



i(p,q) , i{q,p) Hp,q) 



40 



20 



f 

18 



(55) 



d(p,q)e- tx ^\ 

e(p + l,q) +e(p,q- 1) 



■e(p-l,g+l), 



(2p + q)(p + 2q + 6 p -q + 3)/9. 



e(p, q) 
f(p,q) 

i(p, q) 

7 is the eigenvalue of the cubic Casimir operator of 
SU(3). When these operators are inserted into Hamil- 
tonian terms of the type Tr(U\U2Ulu\), expressions in- 
volving up to eight A's and A^'s appear. The variational 
Hamiltonian matrix is generated by taking matrix ele- 
ments of these expressions in our basis. 

Simplifications occur at weak coupling, which we con- 
sider at this point. An imaportant test of the adequacy 
of our basis is to show that the Kogut-Susskind Hamil- 
tonian assumes the form Hyy at weak coupling. At 
weak coupling, the diffusion time t = g 2 r is small be- 
cause we have seen that r remains finite as g is taken 
to zero. Many terms contribute to the sums for S\ and 
S2 when t is small, and the leading behavior at small t 
is yfsTfSl ~ l/(2Vt) = l/(2g^/¥). Thus, at weak cou- 
pling the non-Abelian terms in Jl.b. can be ignored. We 
use Eq. I|27|l . now with the choice a = r. Then at weak 
coupling, we obtain a representation identical to that in 
Sec. El 

1 d 

Jh,a = Jr,u = —^TFT- ( 56 ) 
ig dR a 

At small t the coefficients in U have the behaviors 



(57) 
(58) 
(59) 
(60) 

that are or- 
1(55)1 can be 



DKo 
2& 



We retain contributions to Tr (iJi U2U3UI 

der t or larger. Therefore the term in Eq. 
ignored. Other contributions are: 



K 
Si 


1 - 


- 2t/3, 


K 




Vt 


WS1S2 


2 ' 


K ~ 




25t 2 


Si. 




12 ' 


K 3 


t 


2S 2 


~ ~8 



Tr 



'Ko 
Sx 



3-8t, 



(61) 



Tr 



KoV ( -iK x ~ 

Si j 



AbA c 



(X&+4J 



x (A 2 , c + 4j =-±(2i? 1 i? 2 ), 



Tr 



f) 



4.6^1, 



t 9 t r 

-47^ + 4 



2t. 



(A 1>b f + (Al b f 

Using these expressions and t/r = g 2 , we have 

d? 

2 ^ dR a {s,m) 2 



(62) 



+~ ^2 [R a (s,m) + R a (s + m,n) 

s,m<n,a 

-R a (s + ii, m) - R a (s,n)} 2 

+,i ]T{[A a (s,m)] 2 + [4(s,m)] 2 }, 

s,m,a 

where /1 = — 2t. 

Ha takes the form it does because of delicate combi- 
nations of the small t expressions given in Eqs. I|57|) - 
(|5U1) . The coefficient 1/2 of the middle term of Eq. (52}l 
comes from combining Eqs. (|58|l and l|60() . There is no 
additive c-number in the Hamiltonian because there is a 
cancellation between the term —2t/3 in Eq. (|57|l and the 
term — 1/8 in Eq. I|60l) . together with the normal ordering 
constant 8, the number of SU(3) generators. Note that 
it is surprising that the coefficients we find are even ra- 
tional, because they are ratios of sums that are irrational 
in their small-t limits. For example, the small-i form of 
Si is (7r\/3)/(16t 4 ). 

The value /i = — 2r is unreliable because we have not 
included doubly excited color/link operators in Eq. I|54|). 
The reason such operators appear in Eq. I|62|) (with in- 
correct coefficient) is that the Hamiltonian involves prod- 
ucts of the basic operators. For these terms to contribute 
to matrix elements, one or more of the states must be 
"overexcited". For this reason it might be thought that 
it is harmless to use the value fj, = —2t. That is not 
the case, however, because the Hamiltonian possesses a 
ground state only for < \i < 1/ (4r) . 

The most reasonable choice for /1 is /j, = because 
analogous leftover terms like ^4 a (si)A a (s2), Si 7^ s 2 do 
not appear. The coefficient of these terms is completely 
determined by Eq. I|54|) . 

At this point Ha, Eq. Q62[l. is identical to H\y, Eq. 
(|25H . The expressions were derived in completely differ- 
ent ways, however. Hyy was derived by using weak limit 
expressions for the SU(3) operators. When we derived 
Ha, the operators were unmodified. We constructed the 
Hamiltonian to give correct matrix elements on a limited 
basis of 256 states per link, a procedure that can be em- 
ployed at any coupling. It was only when we took g = 
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at the end, that we found Ha = Hw- We also use Hw 
and Ha differently. We simply diagonalized H w to ob- 
tain exact weak coupling results. Ha will be evaluated 
on our basis of states having small N ave and therefore 
minimal contribution from overexcited states for each 
color/link degree of freedom. Our demonstration that the 
weak coupling form of Ha coincides with Hw is therefore 
only a first necessary check on the viability of the use of 
our basis. What remains to be demonstrated is that the 
spectrum of gauge-invariant states we obtain (by extrap- 
olation to g 2 Qa Ve = 0) is a reasonable approximation to 
the gauge-invariant spectrum of Hw- 

At weak coupling, Q 2 is given by Eq. 1)47(1 (with a — > 
r), and in momentum space N is 



ff = ^BWt(k)fi( s )(k). 



(63) 



k,a s— 1 



We can now compute Gb in momentum space. The op- 
erators N and g 2 Q\ provide a restoring term for the 
longitudinal mode, so for Gb we can perform a Bogoli- 
ubov transformation for both longitudinal and transverse 
modes. The parameters of the Bogoliubov transforma- 
tions are 

4 , m(k) _ -,„{ ^V« + W }, 

4»W(k) - ^W + l + ^ j. (64) 

Then 

2 

G B = G B ,o + ^^G^ 1 , T (k)^)t (k)i? ( S ) (k) 

k.a s=l 

+ ^G B , 1 , L (k)4 3 ) t (k)4 s )(k), (65) 

k.a 

where 



G Bfl /L 3 = 16 A /l + 2£rTi(a) (66) 
+8 v /^77T 1 (/3)-12£, 
Gb.i.t = 2v/[l + 2£T][^ 2 (k)+a], 
Gbxl = 2V[^/r][w 2 (k) + /3]. 

The new expressions in these formulas are: 

a = e/8r, (8 = (1 + 2£r)/8»7, (67) 

T 1(2) (7) = ^EK«+7] 1/2( - 1/2) . 



The average quadratic Casimir operator is: 



.9 2 QLe = 4a 



[Ti(/3)-/?T 2 (/?)]• 



(68) 



We also have: 

2 

k,a s=l 

+iV 2 , r (k) [4')(k)£?«(L-k) 
+4 s )t (k)i? ( S )t (L _ k) - 

+ E{7V 1 , i (k)4 3 )t(k)4 3 )(k) 

k,a 

+7V 2 , L (k) [4 3 )(k)£( 3 )(L-k) 
+4 3 )t(k)4 3 )t(L-k)]}, 



(69) 



L 3 VTTW 



I6r , , yTT^T/ ^ 
Ti(a) H -t 2 (a) 



+4, &!(/?) + J^T 2 (/3)-12, 

V £ r V ^ 



2. 'J' 



v 7 ! + 2£t 8iVw 2 (k) +a' 
r^ 2 (k) +a Vl + 2^r 



VI + 2£t 16r v /w 2 (k) +a' 



1 ^ 



1 



^ v v y 4y n v /^ 2 (k) + /3' 



The mean excitation per color/link degree of freedom 



1 £r 



i 



IS 



No 



24L 3 3V1 + 2£t 



2t , , Vl + 2^ ^, 



24t 
1 



(70) 



It can be shown that the probability that a color/link 
degree of freedom is "overexcited" is 

p = 1 2N 2 ve -2J 2 + 3N ave + l ^ 



[(l + iV a „ e ) 2 - J 2 ] 3 / 2 ' 



where 



,/ = =-- 



2t , N VI + 2fr , s , N 

r2 («) (72) 



3V1 + 2£r 



24r 



At weak coupling, fl^.i = g 2 Q 2 — 0. Computation- 
ally, it is convenient to regroup the terms that remain: 



H B = G C - iN', 



(73) 
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where 

2 

G c = [GB,o-^o]+^^[G B ,i,r(k)-eJV l!T (k)] 

k,a s=l 

x^)t(k)^)(k) 

+ ^[G B , 1 ,i(k)-(iV 1 , L (k)] (74) 

k,a 

Gc incorporates the "diagonal" terms of —£N, and has 
the same eigenstates as the generator Gb- The remaining 
term is 

2 

-ZN> =-£]T5> 2 , T (k) [i?i s >(k)i^)(L-k) 

k,a s=l 

+£Wt(k)^)t(L-k)] (75) 
^E%(k) k 3 Hk)4 3) (L-k) 



+4 3 )t(k)i?C 3 )t(L-k) 



At general coupling, there are many terms in Hb that 
are missing in its weak coupling limit. We will discuss 
this case later. 



IV. THE BASIS AND ITS PROPERTIES 

Generator Gc has an infinite number of eigenstates, 
and we can accommodate only a finite number of these in 
our basis. The goal is a basis that allows us to determine 
the spectrum of low energy, gauge-invariant eigenstates 
of Hb , so we include only the low energy spin zero eigen- 
states of Gc- The lowest energy such states are spin zero 
gluon pair states of the form 

|M = A(k)^^)t(i?k)^)t(i?[L-k])|0). (76) 

R,a 

The wave vector here is one of the "seed" vectors of Eq. 
|@5J. The number of such seed vectors is — 1 + (L + 1)(L + 
3)(L + 5)/48. For example, when L = 31, there are 815 
single pair states of each polarization. These states are 
composed of two gluons and thus are charge conjugation 
even. Both gluons have the same spin index so the states 
have even parity. They constitute a basis for computing 
the spectrum of ++ eigenstates. 

The normalization constants depend on the seed mo- 
mentum because as R ranges over the 24 group elements, 
the number of independent states that appear can be 3, 
4, 6, 12 or 24, depending on k. Denote th e number of 
independent states by V] then A(k) — V2V /96. 

When rj is reasonably small (say, 0.1), the single longi- 
tudinal pair states are closely spaced and lie in an isolated 
band. (At 77 = the band width is zero.) There is then 



a gap between these single longitudinal pair states and 
the lowest energy states composed of two longitudinal 
pairs. So, when 77 is small, our basis will consist of the 
single longitudinal pair states and those single transverse 
pair states that have energies within the band described 
above. In general, only a small number of single trans- 
verse pair states falls within the band of single longitudi- 
nal pair states. This is not a serious limitation because in 
lattice gauge theory the objective is always to determine 
the first few gauge invariant states of given spin, charge 
conjugation and parity. 

As 77 is increased, the content of the basis changes. 
Then the spacing of the single longitudinal pair states 
increases, and further single transverse pair states fall 
within the band. Eventually, the highest energy single 
longitudinal pair state has greater energy than the lowest 
two longitudinal pair state. When that happens, we take 
the basis to include only those single pair states (longitu- 
dinal and transverse) that lie below the two longitudinal 
pair threshold. Note that we do not want to increase 77 
too much anyway, because for 77 > 1, r\g 2 Q 2 becomes the 
dominant operator in the generator. 

When we extrapolate to g 2 Q 2 ve — 0, the longitudinal 
polarization pair states in the basis behave quite differ- 
ently than the transverse pair states. The longitudinal 
states all extrapolate to energy zero, where they merge 
with the ground state. Then they become part of the 
continuum of longitudinal states found in the exact weak 
coupling solution. In our basis, these states are promoted 
to finite energy by the additional operators in Ge- 
lt would be pleasant to simply omit the longitudinal 
polarization pairs in our basis. However, the interaction 
term —£N' connects them to the basis vacuum state |0), 
which is included in the basis. The longitudinal pairs and 
transverse pairs became coupled when we introduced the 
operator N to suppress the "overexcited" states intro- 
duced by our use of boson operators. 

The transverse pair states have the same expectation 
of Q 2 as the basis ground state. That is, they are excita- 
tions of energy only. When we extrapolate to g 2 Q 2 ve = 
the transverse states approach finite limits, and they be- 
come gauge invariant: Q 2 \tp) = 0. We will see that these 
energy limits approximate the energies we found for the 
gauge invariant pair states of Sec. [D] 



A. The Ground State Energy Density 

The energy density we study is the expectation value 
£ = (Gc)/L 3 , as a function of g 2 Q 2 ave - We study this 
quantity rather than (Gb}/L 3 because we want to com- 
pare our results with the weak coupling expression, Eq. 
(|48|l . This ground state energy density has no physi- 
cal importance, but an accurate approximation scheme 
should make a reasonable prediction. We have 



£ 



16(1 + gr) 
VTT2fr 



Ti(a) 



(77) 
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FIG. 1: Ground state energy density vs. g 2 Q 2 ve . L = 31, and 
r\ = 0.1. The lower line is the ground state energy density of 
H w , Eq. J3SJ. 



-iv/TT^Fr 2 (a) + J±T 2 ((3), 

t 2 w r\r 

In Fig. ^ we show this ground state energy density when 
L = 31 and rj = 0.1. The parameters £ and r are cho- 
sen so N ave — 0.5, and g 2 Q 2 ve takes a value on the ab- 
scissa of Fig. ^ The lower branch of the cusped curve 
is the ground state energy density produced by our ba- 
sis. To obtain the energy density at the gauge invari- 
ant limit, we fit a quadratic curve to the lower branch 
of the cusped curve, requiring that the curves intersect 
at the points g 2 Q 2 ave =5.02, 5.00 and 4.98. The energy 
density on this curve at g 2 Q 2 ve — is our extrapolant: 
£ = 19.1192, which is to be compared with the true 
ground state energy density £o=19-1008. The fractional 
error is 0.000974. Very similar results are obtained when 
the quadratic curve is fitted to the lower branch of the 
cusped curve at other points. The probability of "overex- 
citation" of a color/link degree of freedom is in the range 
0.11 to 0.12 for the ground states parameters that pro- 
duce the lower branch of the cusped curve. 

Fig. |21 shows the ground state energy density when L 
and rj are larger. The extrapolated and true ground state 
energy densities are £ — 19.2356 and £q = 19.1008 for a 
fractional error 0.00706. 

The ground state energy curve when r/ = .01 resembles 
Fig. 2] but the cusp is more pronounced. In this case, 
the fractional error of the extrapolant is 0.000316 when 
L = 31. We conclude that the results for the ground 
state are stable when rj is scaled by a factor of 50. 



B. The Basis and Excited States 

The numerical data for the excited basis states is most 
easily explained by laying out a typical example. Con- 
sider the basis when L = 31 and rj = 0.1. We take 
t = 0.428886 and £ = 0.0567920, which sets N ave =0.5 



FIG. 2: Ground state energy density vs. g 2 Qa Ve - L = 51, 
and 77 = 0.5. 
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FIG. 3: Excited state energy (above ground state) vs. g 2 Q 2 ve , 
for eigenstates of Gc- Ordinal states 772 to 778 are shown. 



and g 2 Ql ve = 4.000, on the lower branch of Fig. [T] The 
basis consists of 815 single longitudinal pair states, and 
one single transverse pair state. The transverse state en- 
ergy is independent polarization s, so in fact it is doubly 
degenerate, corresponding to the two-fold degeneracy we 
found for the spin zero eigenstates of Hw- These states 
all lie in the energy interval 0.252337 < E < 0.470661 
above the ground state. Our convention is to order the 
states with increasing energy, with the ground state as- 
signed ordinal number zero. Then when g 2 Q 2 ve — 4.000, 
the transverse pair state has ordinal number 777, and its 
energy is 0.447724. 

The basis energies vary with g 2 Q 2 , and Fig. [3] shows 
seven states near g 2 Q 2 ve — 4.000. The longitudinal pair 
states never cross each other, but the transverse pair 
state crosses four longitudinal pair states in the region 
depicted. The level crossings are consistent with the fact 
that the longitudinal pair energies and the transverse pair 
energy are directed toward different limits as g 2 Q 2 ve is 
reduced. 

We next repeat the computations of the energies when 
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TABLE I: Extrapolated first excited tranverse pair energy. 
Exact energy is 0.4047. For each value of g 2 Qa Ve the top 
entry occurs when the transverse pair is near the top of the 
longitudinal pair basis, as indicated by the ordinal number. 



9 Qave- 

V 


4.000 


5.000 


7.000 


0.043 
0.070 
0.095 
0.100 
0.300 
0.500 


0.3805 c 
0.3788 
0.3076 
0.2483 


0.3843'' 
0.3743 
0.3723 
0.3024 
0.2516 


0.3552° 
0.3379 
0.3236 
0.3236 
0.2453 
0.2018 


"Ordinal number 800/816. 
'Ordinal number 804/816. 
c Ordinal number 804/816. 






TABLE II: 
ergy. Exact 


Extrapolated second excited tranverse pair en- 
energy is 0.5723. 


o 2 2 ■ 
V 


4.000 


5.000 


7.000 


0.070 
0.100 
0.115 
0.190 
0.300 
0.500 


0.5501 c 
0.5232 
0.4725 


0.5692'' 
0.5493 
0.5174 
0.4630 


0.5490" 
0.5357 
0.5259 
0.4900 
0.4456 
0.3848 



"Ordinal number 810/817. 
^Ordinal number 793/817. 
c Ordinal number 587/595. 



g 2 Q 2 ave = 4.000± 0.0200, fit quadratic curves through the 
three data points and extrapolate to g 2 Q 2 ve = 0. The 
extrapolated energies of all 815 longitudinal pair states 
are very close to zero. For example, ordinal state 100 has 
energy 0.313099 at g 2 Q 2 ave = 4.000. The extrapolated en- 
ergy is -0.0110660. To merge with the ground state, the 
extrapolated energy would have to be zero. The frac- 
tional error is about -0.035. As the ordinal quantum 
number of the longitudinal pair increases, the extrapo- 
lated energy of the longitudinal pair also increases. The 
range over all 815 longitudinal pairs is from -0.023 to 
+0.017. Thus, to accuracy we might expect, the longitu- 
dinal states merge with the ground state at g 2 Q 2 ve = 0. 
For the transverse state, however, the extrapolated en- 
ergy is 0.3788, which should be compared to the exact 
result 0.4047 (from Hyp) on this lattice. 

The most interesting issue is how close the extrap- 
olated transverse pair energies approach exact results 
when we vary r\ and the g 2 Q 2 ve from which we start the 
extrapolation. These results are given for the two lowest 
energy transverse pair states in Tables HI and UTI Some ta- 
ble entries are empty because when rj is sufficiently small, 
even the lowest transverse pair states lie above the band 
of longitudinal pair states. That is, they are not in the 
basis as we have defined it. The top entry in each column 
is for a value of rj where the transverse pair state falls just 



below the highest energy longitudinal pair states in the 
the basis, as indicated. 

The next question is how to select a prediction, since 
the table entries vary considerably. For a given column 
in the tables we choose the entry with the smallest value 
of rj, since this parameter fixes the strength of the term 
g 2 r\Q\ in Gb, which perturbs the basis states relative to 
the eigenstates of H a- For each column, this selects the 
largest entry in the column. Given that, it is natural 
to select the largest entry in the table as the prediction. 
The fractional error of this prediction for the first excited 
pair state is then (0.3843-0.4047)/0.4047=-0.050. For the 
second excited pair state the fractional error is -0.005, 
and for the third excited pair state the error is +0.002. 

In the next section we will see that the reliability of 
these extrapolations for transverse pair states is crucial 
for our approach to pure SU(3) lattice gauge theory. 

V. EIGENSTATES OF THE HAMILTONIAN 

The basis consists of one pair eigenstates of Gc, to 
which operator we must add the interaction term = —t;N' 
to obtain Hamiltonian Hb- On our basis, the interaction 
term is a rank one operator. It may be written 

3 

-fJV' - -(384)£A(k)]T]riV 2 , s (k) (78) 

s=l k 

x {|k, s )<o| + |o)<M}. 

The sums are over the seed momenta only. 

We can immediately write half the transverse eigen- 
states of Hb- Recall that the states |k, s) with s = 1,2 
are degenerate eigenstates of Gc- This means that the 
orthonormal combinations 

|k,±) = i={|k,l)±|k,2)} (79) 

are also eigenstates with the common eigenvalue. How- 
ever, it is evident that (— £iV')|k, — ) = 0, so the |k, — ) 
states are already eigenstates of Hb- This means the 
predications we extracted from TableslllandllTlare extrap- 
olated eigenvalues of both Gc and Hb for these states. 
The remaining task, then, is to see whether the states |0), 
|k, +) and |k, L), which are connected by operator —£N' 
form gauge invariant eigenstates of Hb that extrapolate 
to the limits in Tables HI and ITU 

Because —£N' is a rank one operator on our basis, it 
is easy to construct eigenstates of Hb despite the large 
size of the basis. Let |n) denote a pair state belonging to 
our basis, and let E n be the corresponding eigenvalue of 
Gc- An eigenstate of Hb will be a linear combination of 
our basis states: 

\tf>) =a|0}+5>». (80) 

n 

Then the eigenvalue equation 

H B \il>) = (Gc-ZN'M (81) 
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leads to the relations 



H(-flv')|o) 

Pn = c s a, 

t — E n 

£a = ^2(0\(-(N')\n)/3 n . 



(82) 
(83) 



Substituting the first expression into the second, we ob- 
tain the eigenvalue equation 



£ — E„ 



(84) 



It is easy to locate the solutions of the eigenvalue equa- 
tion. Suppose the basis eigenvalues E n are labeled so 
they increase monotonically with n. Then as £ increases 
through the interval [E n , E n +i], the right side of Eqn. 
(|84|l decreases monotonically from +oo to — oo. As £ 
varies in this way, the left side increases monotonically. 
Thus, there is one solution of the eigenvalue equation in 
the interval. In addition to these solutions, there is one 
solution below £ — 0, the new ground state, and one 
solution above the largest E n . 

Fig. 21 shows the eigenstates of Hb in the region of 
the spectrum displayed in Fig. |3J again on a lattice with 
L = 31, T) = 0.1. In keeping with the quantum me- 
chanical result that interacting states do not cross, states 
approach each other closely, but never cross. For this 
reason, there can no longer be a distinct transverse po- 
larization state. However, there is still a trace of the 
transverse basis state. By comparing Figs. |3] and 01 we 
can see how trajectories of eigenstates of Hb bend to 
follow the path of the transverse basis state. Two obser- 
vations explain why this happens. First, when a trans- 
verse basis state crosses a longitudinal basis state, there 
must be an eigenstate of Hb with £ = El = Et- (The 
eigenvector is a linear combination of the two crossing ba- 
sis states.) Between such intersections, eigenvalue equa- 
tion l|84|) still has a solution close the Et because in this 
regime, |(0|(-f A')|T)| 2 |/(0|(-£A')|L)| 2 ~ 0.1. 

The eigenvalues of Hb extrapolate to limits very close 
to zero, and all but two of them behave like nearby 
eigenvalues of Gc- We mentioned that ordinal eigen- 
value 100 of G c extrapolates 0.313099 -> -0.0110660. 
Ordinal eigenvalue 100 of Hb extrapolates 0.313132 — > 
-0.0110627. 

However, two eigenvalues are bizarre. At g 2 = 4.000, 
the new ground state is £q = —64.9678. Note that the 
interaction —£N' is a traceless operator on our basis, so 
the sum of eigenvalues £ is the same as the sum of the 
E n . This suggests that a very positive eigenvalue is to be 
found. Sure enough, £^w = +65.3527. The next lower 
eigenvalue is £$15 = 0.47064, sandwiched between _E 8 i5 
and E$ 16 . 

The very positive eigenvalue is an artifact of the cutoff 
of the basis. For example, if we omit basis state 816, then 
it is £gi5 that is very large. Because of the tracelessness 
of the interaction, the cutoff of the basis also induces the 
very negative £q. 




FIG. 4: Excited state energy (above ground state) vs. g 2 Qa Ve 
for eigenstates of Hb Ordinal states 772 to 777 are shown. 



This pathology largely disappears when we extrapo- 
late. For £ , the extrapolation is —64.9678 — * 0.134561; 
for £si6i the extrapolation is +65.3527 — * —0.137485. 
The fractional deviation of the extrapolated ground state 
energy from zero is +0.002. We will take zero to be the 
extrapolated energy for all these energy levels. 

But if all energy levels extrapolate to zero within ex- 
pected accuracy, what became of the gauge invariant ba- 
sis state that extrapolated near energy 0.4047? At this 
point, the evidence we have of its presence is that eigen- 
states of Hb take turns following the trajectory of the 
transverse basis state. 

This picture is clarified when we follow the "locus of 
gauge invariance," as we vary g 2 Q 2 ve - The expectation 
value of the operator g 2 Q 2 in the eigenstates of Gc or 
Hb may be expressed 



(IsV|) = W m . + a 



(85) 



A gauge invariant state results when we take the limit 
g 2 Qave —> and find A = in the limit. 

Consider A for the basis states, the eigenstates of Gc- 
A transverse basis state has A = at all values of q 2 Q 2 ve , 
and so when we extrapolate this variable to zero, the 
transverse basis states become gauge invariant. The basis 
states |k, — ), share this property, and we have seen that 
they are eigenstates of Hb as well as Gc- 

Among the other eigenstates of Hb, however, no one 
state has A = 0, at nonzero g 2 Q 2 ve For these states, 
A varies with g 2 Q 2 ve - We find there is an abrupt dip 
in A where the state follows the transverse basis state, 
that is, when g 2 Q 2 ve is adjusted so £ ~ Et- This dip is 
explained by Eq. i|82[l . which shows that the amplitude 
of the transverse state in the eigenstate is particularly 
large in this region of g 2 Q 2 ve - 

This behavior is illustrated in Fig. [SJ A dips to a min- 
imum less than 2% of its asymptotic value. Fig. 0] con- 
firmes that the dip occurs when g 2 Q 2 ve is in the interval 
where ordinal state 773 follows the transverse basis state 
trajectory. 
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FIG. 5: A vs. g 2 Ql ve . This is for ordinal state 773. The dip 
occurs where the state follows the trajectory of the transverse 
basis state. 
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FIG. 6: A vs. g 2 Qave- Dips for ordinal states 772 (on right) 
through 777 (on left) are superposed. 



As successive states follow the transverse basis state 
trajectory, the "locus of gauge invariance" passes from 
one eigenstate to the next, always close to the transverse 
basis state. This transition process is shown in Fig. El 
where the dips for states 772 through 777 are shown to- 
gether. Note that dips in A intersect where A is half 
its asymptotic value. At an intersection, the normal- 
ized probability of finding the transverse basis state in 
the eigenstate, |/?t| 2 , is about 1/2 for each of the "dip- 
ping" states. (The sum of the probabilities over all the 
eigenstates of Hb is one.) At an intersection of dips, the 
"locus" passes from one eigenstate to the next. 

Thus, the "locus of gauge invariance" follows the trans- 
verse basis state. This conclusion is guaranteed by Eq. 
(|82|l . The extrapolation of this picture to g 2 Q 2 ve = 
is provided by Table [I] This is as close as we can come 
to the verification that there is a second spin gauge in- 
variant state at energy 0.3843, thereby confirming double 
degeneracy of the spin zero gauge invariant states. 



VI. CONCLUSION 

We now consider how the methods we have developed 
can be applied to pure SU(3) lattice gauge theory at gen- 
eral coupling. In this case, all the terms in Eqs. I|53|) 
and (|54|) contribute to Ha- The "chromoelectric" oper- 
ators J 2 each produce terms which are a product of up 
to four A's and A^'s. The "chromomagnetic" operators 
TrlUllWW) each produce terms which are a product of 
up to eight. The operator g 2 Q 2 produces terms which 
are a product of up to four such factors. 

Terms having more than two factors are interaction- 
type, and are not present in the weak coupling limit we 
have been studying. 

The generator Gb is now constructed from the terms 
in Ha and g 2 Q 2 that are no more than quadratic in A's 
and A''s, using Eq. i|19fl . The basis is chosen from the 
low energy, one pair eigenstates of this generator. 

The main complication beyond those we have studied 
comes from the fact that the expression for Hb, Eq. 12UII . 
now has interaction terms on the right side. When these 
terms are added to —£,N, the resulting interaction term is 
no longer a rank one operator. Therefore, the construc- 
tion of the eigenstates of Hb is more complex than it is 
for the weak coupling case we have studied. 

Interaction terms may pose a technical challenge, but 
they are required if our approach is to produce states like 
glueballs on the lattice. To see this, note that Eq. i|13|) 
implies that a pair state is an excitation in which the 
amplitude for the excitation of two links has the spatial 
dependence e 2,I ' ik As / i j where As is the distance between 
the links. The magnitude of this amplitude is indepen- 
dent of |As|. In a glueball, on the other hand, the am- 
plitude falls off in magnitude when | As| is more than the 
diameter of the glueball. Such behavior can be approxi- 
mated when the eigenstates of Hb are superpositions of 
many pair states, which in turn requires that there be 
interaction terms between pair states. 

One concern is that there may not be enough basis pair 
states to superpose to make a glueball. There are always 
many longitudinal polarization pair states in a basis, but 
that is not true for transverse pair states. An extreme 
example is provided by the case we have used extensively 
in this paper, a lattice with r\ = 0.1 and g 2 Q 2 ve ~ 4.0. 
For these parameters, there are 815 longitudinal polariza- 
tion pair states and just one transverse pair momentum. 
(There are two degenerate states, one for each transverse 
polarization.) To get more transverse states in the ba- 
sis, we increase 77, which increases the energy interval 
between longitudinal pair states. Table ITTll lists the num- 
ber of longitudinal and transverse pair states for different 
values of g 2 Q 2 ve , when 77 = 0.5. These results are for the 
weak coupling generator, but they suggest that it will be 
possible construct a basis with many transverse polariza- 
tion states. 

It should be remembered that adding states to a basis 
can only improve the accuracy of the approximate eigen- 
values we compute. Thus we are not obliged to slavishly 
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TABLE III: Number of longitudinal (L) and transverse (T) 
seed momenta in bases having different g 2 Qa Ve . r\ = 0.5. 



9 Qave- 


4.000 


6.000 


8.000 


10.000 


L 


108 


113 


113 


112 


T 


4 


9 


18 


30 



follow the rules we have set forth for a basis. One can 
add further transverse pairs having energies greater than 
the highest energy longitudinal pairs to improve the "lo- 
calization " of glucballs. One of the attractive features 
of a variational basis is that ideas like this can be tried. 

These considerations set the stage for the exploration 
of the general coupling case, which is the goal of this line 
of research. However, the exploration lies beyond the 
scope of this paper. 
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